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Abstract 

We observe 441 Kl — > e'^e~e'^e~ candidate events with a background of 4.2 events and 

measure B{Kl e+e'e+e- ) = (3.72±0.18(stat) ±0.23(syst)) x lO"* in the KTeV/E799II 

experiment at Fermilab. Using the distribution of the angle between the planes of the 

e+e~ pairs, we measure the CP parameters /3cp = —0.23 ± 0.09(stat) ± 0.02(syst) and 

7cp = —0.09 ± 0.09(stat) ± 0.02(syst). We also present the first detailed study of the e+e" 

invariant mass spectrum in this decay mode. 

PACS numbers: 13.20.Eb, 14.40.Aq, 13.40.Hq, ll.30.Er 
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The rare decay Kl e^e~e~^e~ proceeds via a two-virtual-photon intermediate state 
with internal photon conversions to e^e~ pairs. This permits a measurement of the the 
-^l7*7* form factor, which is necessary to understand the long distance contributions to 
other rare Kl decays, in particular Kl /^^A*^ IQ' whose short distance processes are 
sensitive to the CKM matrix element Vtd- The QED prediction |Q, neglecting radiative 
corrections and a form factor, for the ratio B{Kl e~ e^) /B{K l — > 77), together with 
the experimental measurement of B{Kl — >■ 77) 0, gives a branching ratio of (3.65 ±0.09) x 
10~^. The most precise previous measurements with Kl e^e~e^e~ were based on 27 
candidate events 0|. We report here a measurement of the branching ratio of the rare 
decay Kl — > e^e~e^e'^ at the Fermilab experiment E799II. We also present a study of CP 
symmetry in this decay and the first analysis of a KL'y*'y* form factor. 

The data were collected in 1997 with the KTeV detector at Fermilab. The detector com- 
ponents are described in . During E799II running conditions, 800 GeV protons struck a 
BeO target; collimators and magnets defined two nearly parallel Kl beams, which decayed 
in a 65 m long region held under vacuum. The decay region was surrounded by photon veto 
detectors designed to detect particles escaping the fiducial volume of the detector. A charged 
particle spectrometer, consisting of four drift chambers and an analyzing magnet, was used 
to determine the charge, momentum, and trajectories of particles. A pure Csl calorimeter 
with 3100 crystals was used for photon detection and charged particle identification. 

The four-track trigger used to collect this dataset is described in detail in 0. In the 
offline analysis, we require two positively charged tracks and two negatively charged tracks 
from a common vertex. At least three tracks must strike the Csl; the fourth may strike the 
Csl or pass through one of the Csl beamholes. A track is identified as an if its E/P, the 
ratio of the energy measured by the calorimeter to the momentum measured by the charged 
spectrometer, is between 0.9 and 1.1. We require the energy of the clusters deposited in the 
calorimeter to be greater than 2 GeV. We also place cuts on the track and vertex quality and 
require the photon veto energies to be no more than the level of typical accidental activity. 

To select Kl e~^e~e~^e~ events we require the decay vertex to be within the fiducial 
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region of 95 — 155 m from the target and the total momentum of the kaon in the lab to 
be between 25 — 215 GeV/c. Figure |1| shows the distribution of the P^^, the square of 
the component of the total momentum of the daughter particles (e+e^e+e^) transverse to 
the kaon line of flight, versus the four body invariant mass for the data sample and for a 
Monte Carlo simulation. We define the signal region by Pt^ less than 300 (MeV/c)^ and 
|Me+e-e+e- " Mko\ Icss than 30 MeV/c^, with 90% efficiency. The liberal cut on Mg+g-e+e- 
was chosen to retain radiative events Kl e~^e~e~^e~{j) with an undetected soft photon. 

A source of background to Kl e^e~e^e~ is Kl — > e+e~7 or Kl — > 77, with one 
or two photon conversions in the detector material. The photon conversion probability in 
the material upstream of the first drift chamber in our detector is (2.74 ± 0.11) x 10^^ [0. 
Removing events that have a minimum track separation at the most upstream drift chamber 
of less than 1 mm reduces this background to 3.7 ± 0.3 events, estimated from Monte Carlo. 
Another source of background is Kl Tc^e^u e^e~ where the pion misses the calorimeter 
because of the beamholes and is incorrectly assumed to be an electron. The kinematic 
limit for the invariant mass with the pion misidentified as an electron is 478 MeV/c^. Since 
B{Kl vr^e^z/e^e") is not measured, we estimate it from the branching ratio of the parent 
decay Kl — ^ ir^e^u 7. This background is seen in the lower mass region in Fig. |^ (top) 
and is estimated to be 0.5 ± 0.5 events. After all selection criteria are applied, we observe a 
total of 441 Kl e+e^e^e~ candidate events and estimate the total background to be 4.2 
events, as described above. The distribution of the invariant mass of e'^e^e'^e^ for the data 
and the Monte Carlo is shown in Fig. 

The distribution of the e^e~ invariant masses in the Kl e^e~e~^e~ decay reveals the 
internal structure of the Kl"^*"^* vertex. A related decay, Kl e~^e~j, probes the Kl'J'J* 
vertex; the form factor for that decay has been parametrized by Bergstrom, Masso, and 
Singer (BMS) [|,|, 
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where x = M^e^ /Mk^ and Mk, Mp, Mk*, and are the invariant masses of the 
corresponding mesons. The parameter ax* describes the relative strength of an intermediate 
pseudoscalar decay amphtude and a vector meson decay amphtude. The constant C is 
determined using various couphng constants and is equal to 2.3 [1,0. The BMS model 
predicts \aK*\ = 0.2 — 0.3; note that ax* ~ 0.3, rather than zero, best approximates a 
pointlike form factor {f{x) = 1) for this analysis. Since there are two e~^e~ pairs in this 
decay, we use a factorized expression for the form factor. 

We fit for the BMS parameter ax* from our 441 event sample using a likelihood 
function as in [|T^,^] with one parameter. The likelihood function is based on the 
Kl —>■ e^e~e^e~ matrix element which is a function of all kinematic variables. Our 
Kl —>■ e^e^e^e^ Monte Carlo simulation uses the matrix element calculation found in 
and does not include radiative corrections. Since radiative corrections could potentially 
have a large effect on the form factor measurement, in reality we measure an effective pa- 
rameter, a^* which takes into account both the radiative effects and the form factor. Using 
the likelihood function we obtain aff, = — 0.14 ± 0.16(stat) ±0.15(syst). Figure |^ (a) shows 
the X distribution for data and for Monte Carlo using our measured form factor and using 
a pointlike form factor. The bin-by-bin ratios of data/Monte Carlo (Fig. ^ (b) and (c)) 
further demonstrate the difference between our measured value of aff* and f{x) = 1. The 
systematic uncertainty of 0.15 is due to the detector acceptance, and is dominated by the 
uncertainty in reconstruction of the Z position of the vertex and the calorimeter energy. Our 
result is in agreement with the recent measurement of = — 0.15 ±0.06(stat) ±0.02(syst) 



for the decay e^e 7 [|T^, also ignoring radiative corrections ^. Note that the effect 

of radiative corrections on ax* could well be different in the two modes. 
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We also study the form factor parametrization of D'Ambrosio, Isidori and Portoles (DIP) 



T^j . This parametrization of the ^^^^7*7* form factor is a function of two real coefficients of 
linear and quadratic terms (odip and Pbip, respectively) and is relatively model independent. 
Using a factorization model [jl2[ the expectations for odip and /3dip from a p form factor 
are —1.2 and 1.4, respectively. Due to the dominance of low Meg, e^e^e+e^ is much 

less sensitive to the second order term (/?dip) than the decay Kl e+e^/i+/i^ for a similar 
sample size; thus, only the first order term (aoip) is relevant in this analysis. We measure 
a^ip = —1.1 ± 0.6(stat). This value is in agreement with our measured value of since 
the relation between the two is simply odip ~ — 1 + 2.8aK* (obtained by a Taylor expansion 
of the BMS form factor). 

A first order Taylor expansion of the Kl — ^ e^e~e^e~ form factors yields a generic linear 
form factor (1 + aTayior(a;i + 2:2)), with axayior ~ 0.42 — 1.2ai^. ^ — aDip/2.4. We measure 
an effective parameter alayior = 0-5 ± 0.3(stat). 

We measure the branching ratio of —>■ e+e~e+e~ by normalizing to a sample of 49 
809 Kl n^TT^TT^ decays (vr^ refers to the Dalitz decay 7r° e'^e~7). These events 
were collected using the same trigger as for Kl e~^e~e^e~ . Offline selection criteria 
are similar to those used in the Kl — > e~^e~e^e~ analysis. We require four photons, the 
eight-body (e+e~e+e~7777) Pt^ to be less than 800 (MeV/c)^ and lMe+e-e+e-7777 — Mko\ 
to be less than 20 MeV/c^. The minimum track separation at the most upstream drift 
chamber is required to be greater than 1 mm to reduce conversion background, such as 
Kl 71^71^71^ with a photon from one of the non-Dalitz vr^'s converting to e^e~. The 
background to Kl — > tc^tc^tc^ from photon conversions in detector material is estimated to 
be 844 ± 66 events. Another background is Kl — > ti^ti^tc^l) ('^dd refers to the double Dalitz 
decay —>■ e+e^e+e^ ) since both have the same eight-particle final state. There are 24 
possible combinations of e^'s and 7's that yield Kl tt^ti^tt^ and 3 combinations that yield 
Kl n^n^n^L) ■ We separate the two modes based on the best for the corresponding 
vr'^ decay hypothesis. In addition, we require M^^ to be between 127.5 — 142.5 MeV/c^, 
Me+e-7 between 127.5 — 142.5 MeV/c^ and we exclude the region of Me+e-e+e- between 



127.5 — 142.5 MeV/c^. The number of Kl t['^ti^ti^j~)jd background events is estimated to 
be 21 ±3 events. With these selection criteria, we estimate a total of 865 background events 
for the normalization mode Kl —>■ tt^tt^tt^ . 

The Kl — > e'^e~e'^e~ Monte Carlo, used to estimate detector acceptance, incorporates 
the BMS form factor with a^, = —0.14, our measured value. For the normalization mode 
Kl — s> Tr^TT^jTT^, , a QED calculation of radiative corrections exists in literature [|l^ and is 
included in our simulation. We measure the branching ratio of Kl e~^e~e~^e~ to be (3.72± 
0.18(stat)±0.10(syst)±0.20(norm)) x 10~^ ("norm" refers to the measurement uncertainty in 
the normalization branching ratio 0), consistent with the theoretical prediction 3.65 x 10^® 
0] and the most precise published experimental result (3.96±0.78(stat) ±0.32(syst)) x 10^^ 
1^. Varying a^* by our measured uncertainty of ±0.23 yields a systematic uncertainty in 
the Kl — ^ e^e^e^e^ branching ratio of 0.9%. The addition of final state bremsstrahlung 
in the Monte Carlo simulation does not affect this result significantly. We assign a 1.1% 
error due to limited Monte Carlo statistics and a 2.2% error due to uncertainties in detector 
acceptance, which is dominated by the uncertainty in the track reconstruction efficiency. 
The 5.5% uncertainty in the Kl — > tt^tt^tt^ branching ratio |^ results in an additional 
systematic error. The combined systematic error is 6.1%. Our measurement of the ratio 
B{Kl e+e-e+e-)/B(7rO ^ 6+6-7)^ = (2.59 ± 0.12(stat) ± 0.07(syst)) x 10"^ is unaffected 
by the large uncertainty in the vr^ branching ratio measurement. 

We search for CP violation in Kl — > e~^e~e~^e~ by studying the angle between the 
planes of the two e+e^ pairs in the kaon rest frame. We fit the distribution to |^ 
dTiKL -> e+e-e+e-) ^ ^ ^ ^^^^20) + jcp sin(20). 
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where e is the CP violating mixing parameter and r (which is approximately unity) is 
the ratio of the amplitudes of Ki and K2 decaying to e"'"e~e"'"e~. Ignoring CP violation. 



i.e. neglecting e, Eq. IT] reduces to the Kroll-Wada formula WM, where the constant B 



is — (+)0.20 for an odd(even) CP eigenstate. With the introduction of CP violation the 



coefficient j3cp of the cos(20) term differs from B and tlie sin(20) term, 7cp, no longer 
vanislies. Unfortunately, there are no predictions for the constant C, and therefore for 7cp. 

We measure the CP parameters [3cp and 7cp using the distribution of the angle 0. For 
the purposes of this measurement, we choose the e^e~ pairing that minimizes the product 
of the two e^e~ invariant masses. In reality, both pairings contribute to the matrix element 
since there are identical particles in the final state; the pairing we choose is the dominant 
contribution. To optimize simultaneously the detector resolution of and the statistical 
error, we require the invariant masses of the e^e~ pairs, Meei and Mee2, to be greater 
than 8 MeV/c^, retaining 264 events. Our expectation for B becomes —0.25 as a result of 
discarding low Mge events. We measure I3cp = —0.23 ± 0.09(stat) ± 0.02(syst) and 7cp = 
—0.09 ± 0.09(stat) ± 0.02(syst) after correcting the distribution for detector acceptance 
(Fig. H). We estimate a systematic error of 0.014 and 0.021 in [3cp and 7(7P, respectively, 
due to acceptance uncertainties and an error of 0.008 in both (3cp and 7cp due to detector 
resolution. We also set a 90% confidence level (CL) limit of |7cp| < 0.21. These results are 
consistent with the hypothesis that the decay proceeds predominantly through the CP = —1 
{K2) state. These parameters are insensitive to the effective form factor and to the inclusion 
of final state bremsstrahlung in the Monte Carlo. 

In summary, we observe 441 Kl — ^ e^e~e^e~ candidate events with an estimated back- 
ground of 4.2 events. We measure B{Kl -> e+e-e+e" )= (3.72 ± 0.18(stat) ± 0.23(syst)) x 
10~^. We also study the form factor for the first time in — * e^e^e+e^ . Our mea- 
surements of af.(BMS) = -0.14 ± 0.16(stat) ± 0.15(syst), a%> = -1.1 ± 0.6(stat) and 
Q^ilyior — 0-5 i 0.3(stat) describe the combined effects of a form factor and of radia- 
tive effects. We measure the CP parameters Pcp = —0.23 ± 0.09(stat) ± 0.02(syst) and 
7pp = -0.09 ± 0.09(stat) ± 0.02(syst) (Mge > 8 MeV/c^), thus limiting |7cp| to less than 
0.21 at 90% CL. 
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FIG. 1. Top: The distribution of Pt vs. Mg+g-g+g- for candidate Kj^ — ^ e^e e^e events 
after all cuts except the Pt^ and Mg+g-g+g- cuts. There are 441 events in the signal region defined 



by the box. Bottom: The distribution of Pt vs. Mg+g-g+g- for reconstructed Monte Carlo 
simulated events, scaled to the data statistics. The filled circles represent the signal Monte Carlo 
and the open circles represent the — > ir^e^v e^e~ Monte Carlo. The box defines the signal 
region with an efficiency of 90%. 
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FIG. 2. The M 



e"'"e e'^e 



distribution after all cuts except the invariant mass cut. The dots 



represent the data and the histogram represents the Monte Carlo simulation. The arrows indicate 
an intentionally wide mass window chosen to retain the low-side radiative tail visible in the data. 

The decay Kl — > -K^e^v e+e" is seen in the lower mass region. 
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FIG. 3. (a) The x distribution for data and for Monte Carlo using our measured value of a^* 
and using f{x) = 1. Since there are two e'^e~ pairs in this decay, this is a double entry plot. The 
data/Monte Carlo ratio is shown for our measured value of a|^* (b) and for a pointlike form factor 
(c). We measure a^^* = —0.14 it 0.16(stat) it 0.15(syst). Recall that ax* ~ 0.3 best approximates 
a pointlike form factor. 
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FIG. 4. The distribution of the angle (j) between the planes of the two e^e pairs. In the bottom 
figure, the (f) distribution from 90° — 180° is folded into 0° — 90°. We fit these distributions to Eq. H 
and measure ficp = -0.23 ± 0.09(stat) ± 0.02(syst) and 7cp = -0.09 ± 0.09(stat) ± 0.02(syst). 
The dashed line shows the CP = +l{Ki) prediction and further confirms that the decay proceeds 
predominantly through the CP = —1 {K2) state. 
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